Local three-nucleon interaction from chiral effective field theory 
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The three-nucleon (NNN) interaction derived within the chiral effective field theory at the next- 
to-next-to-leading order (N 2 LO) is regulated with a function depending on the magnitude of the 
momentum transfer. The regulated NNN interaction is then local in the coordinate space, which 
is advantages for some many-body techniques. Matrix elements of the local chiral NNN interaction 
are evaluated in a three-nucleon basis. Using the ab initio no-core shell model (NCSM) the NNN 
matrix elements are employed in 3 H and 4 He bound-state calculations. 

PACS numbers: 21.60.Cs, 21.30.-x, 21.30.Fe 



I. INTRODUCTION 



Interactions among nucleons are governed by quantum chromodynamics (QCD). In the low-energy regime relevant 
to nuclear structure, QCD is non-perturbative, and, therefore, hard to solve. Thus, theory has been forced to resort 
to models for the interaction, which have limited physical basis. New theoretical developments, however, allow us 
Q ■ connect QCD with low-energy nuclear physics. The chiral effective field theory (xEFT) [l[ provides a promising 
bridge. Beginning with the pionic or the nucleon-pion system Q one works consistently with systems of increasing 
nucleoli number @, 0j [3 ■ One makes use of spontaneous breaking of chiral symmetry to systematically expand the 
strong interaction in terms of a generic small momentum and takes the explicit breaking of chiral symmetry into 
account by expanding in the pion mass. Thereby, the NN interaction, the NNN interaction and also 7rN scattering are 
related to each other. At the same time, the pion mass dependence of the interaction is known, which will enable a 
connection to lattice QCD calculations in the future Q. Nuclear interactions are non-perturbative, because diagrams 
with purely nucleonic intermediate states are enhanced Therefore, the chiral perturbation expansion is performed 
j£ | for the potential (note, however, the discussion in Refs. @, H, that points out some potential inconsitencies of 
. ■ this approach). Solving the Schrodinger equation for this potential then automatically sums diagrams with purely 
nucleonic intermediate states to all orders. The xEFT predicts, along with the NN interaction at the leading order, an 
NNN interaction at the 3rd order (next-to-next-to- leading order or N 2 LO) [J [Hi EH j and even an NNNN interaction 
at the 4th order (N 3 LO) [l2l |. The details of QCD dynamics are contained in parameters, low-energy constants 
(LECs), not fixed by the symmetry. These parameters can be constrained by experiment. At present, high-quality 
NN potentials have been determined at order N 3 LO [Hj]. A crucial feature of %EFT is the consistency between the 
NN, NNN and NNNN parts. As a consequence, at N 2 LO and N 3 LO, except for two LECs, assigned to two NNN 
jjj ■ diagrams, the potential is fully constrained by the parameters defining the NN interaction. 

It is of great interest and also a challenge to apply the chiral interactions in nuclear structure and nuclear reaction 
calculations. In a recent work [l4|, the presently available NN potential at N 3 LO [l3[ and the NNN interaction at 
N 2 LO 0, HH have been applied to the calculation of various properties of s- and p-shell nuclei, using the ab initio 
no-core shell model (NCSM) [H, HBl, up to now the only approach able to handle the nonlocal %EFT NN potentials 
for systems beyond A = 4. In that study, a preferred choice of the two NNN LECs, cd and ce, was found and the 
fundamental importance of the %EFT NNN interaction was demonstrated for reproducing the structure of mid-p-shell 
nuclei. In a subsequent study, the same Hamiltonian was used to calculate microscopically the photo-absorption cross 
section of 4 He fl7j |. 

The approach of Ref. [3 differs in two aspects from the first NCSM application of the xEFT NN+NNN interactions 
in Ref. [l8[, which presents a detailed investigation of 7 Li. First, a regulator depending on the momentum transfer in 
the NNN terms was introduced which results in a local xEFT NNN interaction. Second, the 4 Hc binding energy was 
not used exclusively as the second constraint on the cd and ce LECs. 

A local NNN interaction is advantages for some few- and many-body approaches because it is simpler to use. At 
the same time, it is known that details of the NNN interaction are important for nuclear structure applications. 
For example, the Urbana IX and the Tucson-Melbourne [2(| Hi], [22| NNN interactions perform differently in 
mid-p-shell nuclei [l6|, HH, [24| although their differences appear to be minor. In the Green's function Monte Carlo 
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(GFMC) calculations with the AV18 NN potential [25[, the best results for p-shell nuclei up to A = 10 are found using 
the Illinois NNN interaction that augments the Urbana IX by a two-pion term from the Tucson-Melbourne NNN 
interaction and by three-pion terms that in the xEFT appear beyond the N 3 LO [26|, [27} ■ Contrary to the Illinois 
NNN interaction, the xEFT NNN interaction features the above mentioned consistency with the accompanying NN 
interaction. Still, interestingly, we found that the nonlocal xEFT NNN interaction used in Ref. [l8[ and the local 
XEFT NNN interaction employed in Ref. [l4| differ to some extent in their description of mid-p-shell nuclei with the 
latter giving results in a better agreement with experiment. Therefore, it is important to pay attention to the details 
of the NNN interaction and test different possibilities. 

It is the purpose of this paper to elaborate on the details of the local xEFT NNN interaction used in Refs. [3, IT?} 
and present its matrix elements in the three- nucleon basis. Technical details of dealing with NNN interactions were 
investigated in many papers [H, [29|, [3(| [U, [H, [H, [34| . A new feature in the present work is the use of xEFT contact 
interactions and a focus on the application within the ab initio NCSM. In particular, we demonstrate the binding- 
energy convergence of the three-nucleon and four-nucleon systems with the xEFT NN+NNN interactions using the 
ab initio NCSM. In Sect. HH the local xEFT NNN interaction is discussed and compared to the nonlocal version of 
Ref. [HI]. Its three-nucleon matrix elements are given term by term. In Sect. IIII| the 3 H and 4 He binding energy and 
radius calculation results using the N 3 LO %EFT NN interaction of Ref. [l3| and the local xEFT NNN interaction arc 
given. Conclusions are drawn in Sect. IIVI 



II. LOCAL X EFT NNN INTERACTION AT N 2 LO 

The NNN interaction appearing at the third order (N 2 LO) of the xEFT comprises of three parts: (i) The two-pion 
exchange, (ii) the one-pion exchange plus contact and the three-nucleon contact. In this section, we discuss all the 
parts in detail and present the three-nucleon matrix elements of all the terms. For the two parts that contain the 
contact interactions, we also discuss in detail the impact of different rcgularization schemes. 

A. Three-nucleon coordinates 

We use the following definitions of the Jacobi coordinates 

ei = ^fa-rd, (i) 

6 = J\ (\(n+r2)-r 3 ) , (2) 



and associated momenta 



7?i = -/=(pi -V2) , (3) 

^2 = \/| ( 7T(Pl+P2) -P3j • (4) 



3 V2 

We also define the momenta transferred by nucleon 2 and nucleon 3: 

Q = = -^|(7fi -7fi) + ^(4 -^2) , (5) 

Q' = p 3 -P3 = -y|(*2 - ^2) , (6) 

where the primed coordinates refer to the initial momentum and the unprimcd to the final momentum of the nucleon. 



B. General structure of three-nucleon interaction and its matrix element 



The NNN interaction is symmetric under permutation of the three nucleon indexes. It can be written as a sum of 
three pieces related by particle permutations: 



W = Wi + W 2 + W 3 



(7) 
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To obtain its matrix element in an antisymmetrized three-nuclcon basis we need to consider just a single term, e.g. 
W\. In this paper, we use the basis of harmonic oscillator (HO) wave functions. However, most of the expressions 
have general validity. Following notation of Ref. [35| , a general matrix element can be written as 

(NiJT\W\N'i'JT) = 3(NiJT|Wi|JVYJT) = 3^2(nlsjt,N 'CJ\\NUT)(n'l' V 'ft' ',Af' C 'J'\\N' Y JT) 

x((nlsjt,Af£J)JT\Wi\(n'l's'j't',Af'£'j')JT) , (8) 

where \NiJT) is an antisymmetrized three-nucleon state with N = 2n + 1 + 2Af + C, i an additional quantum number 
and J and T the total angular momentum and total isospin, respectively. The parity of the state is (— 1)^. The state 
\{nlsjt,NCJ)JT) is a product of the HO wave functions (£i\nl) and (£ 2 |A/"£) associated with the coordinates (UJ 
and ([2|), respectively. This state is antisymmetrized only with respect to the exchange of nucleons 1 and 2, i.e. 
(— l) i+s+ * = —1. The coefficient of fractional parentage (nlsjt,N CJ\\NiJT) is calculated according to Ref. [35|. 



C. N 2 LO three-nucleon interaction contact term 



We start our discussion with the most trivial part of the x^FT N 2 LO NNN interaction, the three-nucleon contact 
term 

WT nt = Ef 2 ■ f 3 6(n - f 2 )S(r 3 - r!) = Ef 2 ■ f 3 ^^L^ J dtfidtfad^d^i^X^I , (9) 

with E = where A x is the chiral symmetry breaking scale of the order of the p meson mass and = 92.4 MeV 

is the weak pion decay constant. The ce is a low-energy constant (LEC) from the chiral Lagrangian of order one. 
The corresponding diagram is shown in Fig. [1] This term was regulated in Ref. [Tl[ by a regulator dependent on the 
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FIG. 1: Contact interaction NNN term of the N 2 LO xEFT. 
sum of Jacobi momenta squared: 

^cont.ENGKMW = Er 2 ■ T 3 ^ J- J d^d^ \H^)F{\ (tt? + 7T 2 ) ; A)F(| « 2 + ^ 2 ); A) <^ 2 I , (10) 
with the regulator function 

F(g 2 ;A)=cxp(- (Z 4 /A 4 ) (11) 

with the limit F(q 2 ; A — > oo) = 1. This was in particular convenient as the calculations were performed in momentum 
space. 

Alternatively, let us consider a regulator dependent on momentum transfer: 

w cont,Q = m .^_±_ 1 |d^d^d 7 rid4|^^)F(g 2 ;A)F(Q' 2 ;A)<^ 1 7r^| 

- Er 2 ■ t 3 J d£d£||i&)2o(>/2&; A)Zo(|^6 + yflU, A)(|*i6| , (12) 
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where we introduced the function 

Z (r;A) = ^J dqq 2 j Q (qr)F(q 2 ;A) . (13) 

This results in an interaction local in coordinate space because of the dependence of the regulator function on 
differences of initial and final Jacobi momenta. An interaction local in coordinate space may be more convenient for 
some methods. In fact, most of the NNN interactions used in few-body calculations, such as the Tucson-Melbourne 
(TM') [H|, HI], Urbana IX (UIX) [H[ or Illinois 2 (IL2) (Hi, are local in coordinate space. 

The two alternatively regulated contact interactions lead to different three-nucleon matrix elements. The interac- 
tion (jTUJ) gives 

((nlsjt,Af£J)JT\W™ nt < ENGKMW \(n'l' S 'j't\Af'£'j')JT) 

i ,,,,,,,,, ^ t+t > + T+i r t f 1 1 r t f i 



2V07T 4 1222J122 1 



x J dn 1 dw 2 w 2 w 2 (-l)^+^R n0 (w 1 , i)i?AAo(7r 2j i)f(!(7r? + 7r2); A) 

x J d7rid^< 2 7 r 2 2 (-l)("'+ AA ')i?, i , K,i)i?^ (<,i)F(iK 2 +^ 2 );A) , (14) 

while the interaction (|12|) results in the following matrix element: 
((nlsjt, NCJ)JT\ w™ ntQ | (n'l's'j't',M'£'J')JT) 

= E65„ti>(r-ir* +T+ *l t 1 { \ U \ i ;)]?jjwH) j -^ j '- j+i+/;+s 

I 2 2 2 J I 2 2 J 

x £(-l)^ 2 {] ]> X s }{j> j X j}{ J c C \ ] (i'Q^0|i0)(£'QX0|£0) 

x / d&d&£$iM6»&)^£(6>&)iW^^ (15) 

In the above expressions, we have introduced the radial HO wave functions R n i with the oscillator parameter b and, 
further, a new function 

Z , x (ri,rr,A) = ^ J dqq 2 jx(qri)jx{qr 2 )F(q 2 ; A) . (16) 

We also introduced the customary abriviation I = y/21 + 1. It should be noted that the two differently regulated 
contact interaction have different tensorial structure. One would perhaps expect that matrix elements of a local 
interaction will be easier to calculate. This is not the case for the discussed contact interaction. From Eq. (fT4|) we 
can see that the term (fTOj) acts only in 5-waves. On the other hand, the local interaction (p~2|) acts in higher partial 
waves as well as seen from Eq. (|15p . We display this schematically in Fig. [2] by breaking the symmetry of the pure 
contact interaction diagram (Fig. [T]) and showing the finite range of the momentum transfer regulated interaction. 
Using F(q 2 ;A — > 00) = 1, Zq(t, A — > 00) = j^?5(r) and Zo } x(ri, f 2 ; A — > 00) = j^8{r\ — r 2 ), it is straightforward 

to verify that in the limit A — > 00 both expressions (fl"4|) and (fl"5|) lead to the same result. For completeness, let 



us note that in Ref. [ll| a matrix element of f\ ■ f-2 was calculated, i.e. (ti • T2 

) = 6<5i f /(-l)*- 1 I I I * j, instead 



of T2 • T3 as we do in Eq. (fT4j) . Either choice lead to identical matrix clement in the three- nuclcon antisymmctrized 
basis ©. This is not the case once we regulate with the momentum transfer. Our choice in (TT2"|) results in the same 
isospin-coordinate structure as that obtained in Ref. [lOl ]. 

D. Transformation of the momentum part of the NNN interaction 

A general NNN interaction term is a product of isospin, spin and momentum parts. In this subsection, we manipulate 
the momentum part. We only consider the case of the regulator function F(q 2 , A) depending on transfcred momentum. 
The momentum part of a general term W± can be schematically written as 

g Kl (\Q\;A)g K2 (\Q'\;A)(Y Kl (Q)Y K2 (Q'))^ , (17) 



■5 
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FIG. 2: Contact interaction NNN term of the N 2 LO xEFT regulated by a function depending on momentum transfer. 



with K\ + Ki even and with Q and Q' defined by Eqs. ((3]) and (j4|), respectively. For coordinates and momenta, Q 
denotes the angular part of the vector Q. A transformation of (fTT]) to coordinate space leads to a local interaction 



(2tt)6 (V3)3 



d^ 1 d^d7fid^| 7 r 1 7r 2 ) 5 K 1 (|g|;A) ffif2 (|Q'|;A)(y Kl (Q)r i , 2 (Q')) (i " ) ('?i4 



■k 1 +k 2 



J dad|* 2 |af 2 }/^ (V2C1; A)/if 2 (i^a + ^fai; (a)*v 2 (^a + vfe)) w • ( 18 ) 



Using ([I]) and ([2]), we note that V%£i = ri—r 2 and -^a + y fa = ri — F3. In the above equation, we have introduced 
a new function using the relation 



which implies 



i K f K {r;A)Y Kk {r) = — / dq ^ r g K (q;A)Y Kk (q) 



f K {r;A) = / dqq 2 j K (qr)gK(q;A) . 



We manipulate Eq. (|18[) first by utilizing the spherical harmonics relation 



K 2 



Yx 2 k 2 (ri +f 2 )= 53 



A' 3 =0 



#3 L 



2/v 2 + 1 
2if 3 



rfvf 2 - lf3 |f 1 +f 2 |-^ 2 (y^(f 1 )y X2 _ if3 (f 2 ))^ 



-A\, 



K 2 



(19) 



(20) 



(21) 



and, second, by the following expansion involving the functions depending on |^a + y § ai : 

/K 2 (|ri+f 2 |;A)|fi+f 2 |-^ =4tt £ /^(n, r 2 ; A)(-l) x Fl Mx {h)Y XMx (r 2 ) , 



with the function Jk 2 ,x(j\, r 2 ; A) given by 



fK 2! x(n,r 2 ; A) = - / dqq 2 j x (qri)jx(qr2) [ drr 2 j (qr)- K 



or, equivalcntly, by 



/K 2 ,x(ri,r 2 ;A) = § duPy (tt) ^^^jp 



/k 2 (y>^+r|-2r 1 r 2 M;A) 



(22) 



(23) 



(24) 
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Using (|21[) and (f2"2"| . the term (JTHJ) is re- written in the form 

K 2 



■K 1 +K 2 



K 3 =0 XYZV 



(2K 2 + 1 
V 2^3 



(7f6) X3 (\/|6)^-^/i f2 ,x(^a^/|6;A) 





y 




{5 


K 





xX 2 K 2 - K 3 K 1 K 2 Y(-l) Kl+Y+z+K 

x(xok 3 o\yo)(xok 2 - /v 3 o|zo)(yox 1 o|yo)(y y (e 1 )yz(6))r ; 

which is convenient for matrix element calculations. 



Y X K 3 
K 2 -K 3 K 2 Z 



(25) 



E. One-pion-exchange plus contact N 2 LO NNN term 

We are now in a position to discuss the one-pion exchange plus contact term that appears at the N 2 LO. Following 
Ref. [ll|, we can write the W\ term contribution as 



Trrl7r_cont n 1 9A -j 



1 



Ml 



S 2 ■ Q'3 3 ■ Q' 



1 



Ml 



o 2 ■ Qa 3 ■ Q 



(26) 



with D = F %R , where ca is a LEC from the chiral Lagrangian of order one. A diagramatic depiction of the second 
term in the parenthesis is presented in Fig. [3l The first term corresponds to the exchange of 2 <-> 3. Using the 
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FIG. 3: One-pion exchage plus contact NNN interaction term of the N 2 LO xEFT. 
regulator dependent on the sum of Jacobi momenta squared of Ref. f]~lj | , this term can be cast in the form 



l7r.cont,ENGKMW 



-D 



1 9A 



(2tt)6 8F; 



— jL- J d^dTfadTfid^l^F^Tr? + 7r|);A) 



XT 2 • T 3 



1 



-B 2 ■ Q'S 3 ■ Q' 



1 



<?2 • Q&3 ■ Q 



_Q ,2 +M% Q 2 + M 2 

On the other hand, with a regulator dependent on momentum transfer, we get 



F{Wi +^);A)«7? a |. (27) 



w 



l7r_cont,Q 



= -D 



1 ffA 1_ 

(2tt)6 8F 2 (^3)3 
1 



XT 2 • T 3 



dTridTradTfidTfa |7fi7? 2 )F 1 (g 2 ; A) 
1 



Q' 2 + M 2 



<?2 • Q'<? 3 • Q' 



M 2 



^2 ' Q(T 3 • Q 



F(Q' 2 ;A){^' 2 \ , 



(28) 



which leads to a term local in coordinate space. We depict the second term in the parenthesis of (|28|) schematically 
in Fig. 21 The first term corresponds to the exchange of 2 <-> 3. This choice of regulation results in spin-isopin- 
coordinate structure that also appears in NNN terms obtained in Ref. [13] ■ We note that a somewhat different 
spin-isopin structure was used for pion-range-short-range NNN terms in Refs. [32[ and [34 [. In Ref. [13] in particular, 
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the a and f operators were associated with the active nucleon 1, i.e. 



w; 



l7r_cont,Q CT1 



-D 



1 9A 1 
(2tt)6 8F 2 (^3)3 

1 



d7fid7r 2 d7r , 1 d7f^|7r 1 7r 2 )F(Q 2 ; A) 
1 



n • r 3 



Q' 2 + M2 



<?1 • Q'(T3 ■ Q' + T\ ■ T2- 



1 • <3<72 • Q 



F(Q' 2 ;A)(^ l7 r 2 | , 



(29) 



This change docs not alter the matrix element of l|27j) in the antisymmetrized three-nucleon basis. It will lead to 
a difference in the matrix element of {28]) ■ However, the dependence on the regulator is a higher order effect than 
the xEFT expansion order used to derive the NNN interaction. Therefore, these differences should have only minor 
overall effect. In fact, we confirmed in nuclear structure calculations such as those described in Ref. [l4[ that impact 
of the choice (|2"T|) or (|2"5|) is small in particular when the natural LECs values are used (\cd\ ~ 1)- However, a more 
significant impact of the choice of the regulator in particular on spin-orbit force sensitive observables is observed in 
the case of the two-pion-exchange terms as discussed in the Introduction. 




FIG. 4: One-pion exchage plus contact NNN interaction term of the N 2 LO xEFT regulated by a function depending on 
momentum transfer. 

Due to the antisymmetry of the three-nucleon wave functions in (|S]), it is sufficient to consider just one term of the 
two in parenthesis in (|27|) and (|28|) and multiply the result by two. Using the first term, the matrix element of ([27)) 
with the regulator dependent on the sum of Jacobi momenta squared is obtained in the form 

l7r_cont,ENGKMW i / / 



({nlsjt,AfCJ)JT\W; 
-- -D^ — ^—-SiqSi'qSsjSs'j'H' 



\{ril's'j't\N'C'J')JT) 



t t' i 

2 2 2 



t i! 1 

1 I T 

2 2 



J-J+T+ 



{ { < 






f 




I 2 2 






J 





K i=0 



2K + 1 
2Ki 



(-1) 



C+K 



| (^—l^n+AT+n' +AT') + (C+C')/2 

( J C i 

K(WW\K0) ij'C'l 

K=0,2 [ I A" 1 

1 XC'iKiOXOlCOXC'OK - KtflXO) 



x 



£ K-K x X 
K 1 £ K 



x J dTridTradTridT^TrjTrlTif ^^(tt!, \)Rmc{^2, ±)F(±(^ 2 + tt 2 2 ); A) 
xi?„, W,i)i?^K,i)A(iW 2 + ^ 2 );Ahf 1 7rf- Kl ^x(7r 2 ,7r 2 ) , 
where we introduced the function 



2 f q 2 ~ K 
9K,x{p,p') = - / dqq 2 drr 2 r2 jo{qr)jx(pr)jx(p'r) 
7T ./ yr + M£ 



which can be alternatively evaluated through 

i f 1 



, \ V 'p 2 +p' 2 -2pp' u 

9K,X(P,P ) = 5 U AuPx ^ ^\^- 2pp 'uf + Ml 



(30) 



(31) 



(32) 



The matrix element (|30p was first derived in Ref. [1JJ 
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For the one-pion exchange plus contact term (|28p with the regulator dependent on momentum transfer, we present 
the matrix element obtained using both terms in the parenthesis of (|28[) . Due to the three- nucleoli wave function 
antisymmetry, both contributions lead to the same result for ((5J). One can take the advantage of this feature and use 
the alternative calculations to check the correctness of the numerical code. First, we take the first part of ([28]) and 
get 

l-7r_cont,Q| a /,//•/./ \f C> T'\ 



((nlsjt,Af£J)JT\Wi* omM \(n'l' V 'ft' ' ,N' '£' ' J')JT) 
_ D |A fT( _ 1)t+t ' + T + i / t t> 1 



l222jl22 J J I 2 2 2 J 



x ^(-l) W2) (1010|ii'0) 53 K 2 K 3 l'£'(l'0K 2 0\l0)(£'0K 3 0\£0) 

K=0,2 K 2 ,K 3 

£ § J 

£' \ J 1 

K 3 1 Z 
1 










f a 2 1 z\ 




J 


5} 


1 1 a 3 a / 



2 ^ i\Ki\^r i\Xv2 tT s nvcwis (WITS ts nvr>\ ts rA J %1 X K 2 

K 3 A K-Ki 



-l) Kl J2{-l) X X 2 {Ki0X0\K 2 0)(K - KrfXOlKsO) 

x 

x J d^ 2 ei^Rni^i,b)R^c(^b)R nl ^^ub)R^ c ,(^,b) 



><(Vki)^(Vi6) X " Xl ^o(v^i;A)/ x , JS -( > /|€i, > /|6;A) , (33) 
with the functions 

/o,x(n,r 2 ;A) = / dqq 2 jx{qri)j X {qr 2 ) q ^77^ , (34) 



^2 J l,./^ g2 + M 2 ' 

and 

1 f 00 f 00 F(k 2 ;A) 

f2.x(ri,r 2 ; A) = — y dqq 2 ix{qri)jx{qr 2 ) J dkk(k 2 - g 2 ) fc2 - ^ , (35) 

which are special cases of (|2"3"j) . An alternative way of evaluating (|55|) is 

-f / AN 1 r 1 j t-> / \ /2(v/ rf+r^-2r 1 r 2 u;A) 

/ 2 ,x(r-i,r 2 ;A) = 5 J_ 1 duifr (u) ^'^.^ , (36) 



with 



Mr;A) = ^/d^^)2^2. (37) 

We note that (|36p is numerically more efficient than (j35"T) . 

Next, we take the second part of (|28[) . which results in a simpler expression for one-pion-exchange plus contact 
N 2 LO three-nucleon matrix element in non-antisymmctrized basis: 

((nlsjt,M£J)JT\Wt™ ont ' Q \(n'l's'j't',M'£'J')JT) 

= -D^-tt'(-iy +t ' +T+ 2 / * n }{ 5 ? i, bj-w^l ; i \ 

^tt L2 2 2JL2 2 J J I222J 

53 a'(-i)W 2) (ioio|a'o)53(- 1 ) v '^(i / o/'o|w)53a 2 (aoa-o|i/o)(ao/:'o|/:o) 



x 

if =0,2 



X 

z 




E^ 2 

z 

J dfid&^fikiKi. b)Rtf C fe, b)R n ' V ^ 1 ,b)R Aric >^ 2 ,b)f K (V2^;A)Z 0!X (^ 1 , ^§6; A) , (38) 
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f K (r; A) = -L / dgg 2 j K (gr) - (39) 



with 

f ~<-^ J ->*WV g 2 +M 2 ' 

which is a special case of ([20]) and ^o.xC^ij^a; A) given by Eq. ([IB)) . Both ([55)1 and are already multiplied by two 
in anticipation of the three-nucleon antisymmetry in the final matrix clement ((SJ . 
For completeness, we also present the matrix element of the second part of (|2!)1) : 

((nlsjt,MCJ)JT\wl % - cont ' Q,rl \(n'l's'j't',WC'J')JT) 
F 2 



x 



| | j J |jj' l 7 l 7'^'f'£'(-l) s+s '+^+ J+£ +- 7+ - 7 ' + 5 

( s s' Jf ] 

^ /^(-l) (/c/2) (1010|^0) < | ± 1 I ^V r (F0Z / 0|Z0)^X 2 (X0iT0|F0)(X0£ / 0|£0) 
:=o,2 I » 1 y x 




f J X J" 

x / d&d6£?£iJM6,&)i^(&>&^ , (40) 

which is still simpler than (|38[) . Again, this matrix element is already multiplied by two in anticipation of the three- 
nucleon antisymmetry in the final matrix element (J5J). The matrix element of the (two-times the) first part of (J^HJ) is 
given by (|33|) multiplied by (— l) t+t + s + s . Due to the three-nucleon wave function antisymmetry, this contribution 
leads to the same result for ([5]) as does which can be taken advantage of in testing the correctness of numerical 
calculations. 

By comparing (|27|) with (|33|) (or cquivalently with ([38]) and also with (|40|) ) we note the different tensorial structure 
of the matrix elements. When the regulator dependent on the sum of Jacobi momenta squared is used, only the / = 0, 
I' = partial waves contribute. This is not the case, when the regulator depending on momentum transfer is utilized. 
At the same time, however, we note that in the limit A — * oo both expressions (|27[) and as well as ([55]) and ([4T)]) 
lead to the same result. 



F. Two-pion exchange N LO NNN terms 

In this subsection, we present matrix elements of two-pion exchange N 2 LO NNN terms. Their schematic depiction is 
shown in Fig. [5] There are three distinct terms associated with three LECs, ci, C3 and C4, from the chiral Lagrangian, 



3 1 2 

FIG. 5: Two-pion exchange NNN interaction term of the N 2 LO xEFT. 



which also appear in the sublcading two-pion exchange in the NN potential. Consequently, values of these LECs 
expected to be of order one are typically fixed at the NN level unlike the case of the previously introduced cd (|26|) 
and ce © LECs whose values needs to be fixed in systems of more than two nucleons. In the present paper, we 
derive only the matrix elements of the two-pion exchange NNN interaction terms regulated by a function depending 
on momentum transfer, i.e. terms that are local in coordinate space. 
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Following Ref . [Ill ] , the W\ part of the c\ term with the momentum transfer regulators can be written as 

w " - fk* ■ f > F ^ A) whrr - fe ■ Q'wrirf <«' 1; A) ■ (41) 

Using results of Subsect. Ill Dt we find for the ci-term matrix element: 
((nlsjt,MCJ)JT\W^- clM \(n'l's'j't',M'C'J')JT) 



f2 n 2 



7T 7T L222JL22 J L222 



(I s j ) { L i J 

y2VR(V0l / 0\l0)(R0£ / 0\C0)Y j (-l) Y Y(Y010\V0) { I' s' f £' | J' 

vr y { V 1 Y J [ R 1 Y 

j 1 



J' J J 



X 

if 3=0 



x (V^i) K3 (V^ 2)1 " K3/l( ^ i;A)/l ^ ( v^ i 'V / ^ 2;A ) ■ ^ 

Here we introduced the functions 

h(r;A) = ^Jd qq ^( q ) q -p^, (43) 

and 

1 f 00 Z" 00 kF(k 2 - A) 

fi,x(ri> r X A )=2^S J o dqq 2 ]x(qri)]x(qr 2 ) J dk fc2 + ^ 2 , (44) 

which are the explicit versions of functions given in Eqs. and (|123"|) . respectively. The function (|4"4")l can be 
alternatively evaluated with the help of the Legendre polynomial: 



/i,x(rx,r 2 ;A) ,I/^d^(«)^^S . (45) 
The Wi part of the two-pion exchange c 3 term is given by [ill ] 

^ " «otIt 2 ■ f3F(<?; A) ^T ' fe ■ ' ^'o^bf" 5 " A) ' (46) 

For its matrix element we find 

((nlsjtMCJ)JT\W^- c3 ' Q \(n'l's'j't\Af'C'j')JT) 

. Afft.,,^ j ( « j j j j « j 1 ) rm , sn _ iy . J+ . + , r./i W 

7T 7T L222JL22 J L222J 

x (-l) (A ' 1+K2)/2 ^i^2(1010|A' 1 0)(1010|X 2 0) 

K 1 K 2 

J2(-l) v+R VR(V0l'0\l0){R0C'0\C0) ^ ^(^OAiOlUO) 



x 

VR Y 



if 2 r 



f 2 ^ + 1 ) " K ^I<3 X 2 (X0K 3 0\Y0)(X0K 2 - K 3 0\R0) 



x 

K 3 =0 



Y X K 3 
K2-K3 K 2 R 



V 2X 3 

x j ^ l ^ 2 ^tR n i^i,b)R NC {^b)R n n'^iMRN'C'{^,b) 

x(x/^a) K nx/^6) K2 " K3 M(V2a;A)/ Jf2 ,x( A /iCi,A/|6;A) , (47) 
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with the function fK 1 {r;A) given by (|3"9"j) . the function fo,x{ri 7 r- 2 ; A) given by (|3~4"|) and the function fe t x (fi, r 2 ; A) 
given by ((35)) . 

Finally, the Wi part of two-pion exchange C4 term is given by [ll| 

^" C4 = ^ • ft >< T3)F(Q 2 ; A)^^. 2 • ■ Q'oi • (Q X (ft^-A—J^; A) , (48) 

and for its matrix element we find 

((nlsjt,M£J)JT\W? w - ci ' Q \(n'l's'j't',N'£'J')JT) 



X 



2 (-l)( Kl+X2 '/ 2 A: 1 A > 2(1010|A' 1 0)(1010|A' 2 0) 



R(V01'0\10){R0£'0\£0) ^ r(Yo.sriO|vo) 



X 

x]Tz 2 (-lf £' I J' 
z \ R 1 Z 





f 




1 z 


1 y| 




j 


?} 


1 A 2 


i? 1 / 



E*?f *' ' 'All I ' 

x 4 [ V Ki Z J [ 1 1 K 4 

K2 





1 


Y \ 




1 




u 


V 


A 4 J 


{i 


A 4 


A} 



X 

if 3=0 



y a a 3 

A 2 -A 3 A 2 R 



<r 3 =o L V 3 / J x I 

x J d^d^lRni(ix^RM^,b)Rnn'^x,b)R^'^2,b) 

^(^^'(^^-^V^Cv^ijA)/^^!,^^) . (49) 



The same functions (|34|) . (|35|) and ([39]l that were introduced in the C3 term enter the C4 term as well. 

We note that the local two-pion-exchange terms appear also in the Tucson-Melbourne NNN interaction [2(J. The 
analogous terms to c\ , C3 and C4 are present in particular in the TM' interaction [2ll . |22| . The TM' parameters are 
denoted by a', b and d with the relation to the above c\, C3 and C4 given by 

a' = , (50) 

b = 4c 3 , (51) 

7T 

d = ^c 4 . (52) 

Further, the regulator function F(q 2 ;A) is chosen in the form 

A 2 - M 2 

^V;A) = ^r- (53) 

This choice allows to evaluate integrals that define the functions Jk analytically. The analytic expressions can be 
found, e.g. in Ref. [29l |. In that paper, the following function is introduced: 

Using the properties of spherical Bessel functions, we can easily find relations between derivatives of Z\ and our fx 
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functions: 



/o(r,A) = - fe(r;A) + H^(r;A)) , 



(55) 



/i(r,A) = -Z[(r;A) (56) 
/ 3 (r,A) = ^'(r;A)-^(r;A) . (57) 

For completeness, we note that a still different notation was used in Ref . [Io| , where functions Ik,i (r; A) were introduced. 
They are related to the Zo(r, A) that we introduced in Eq. (fl~3|) and to the above Zi function (|54|) is as follows: 

Z (r,A) = / 0l o(r;A) , (58) 
Z x (r,A) = / 2 ,o(r,A) . (59) 

In Ref. [36| , the Tucson-Melbourne NNN interaction matrix elements in the HO basis were calculated using a differ- 
ent algorithm than the one used in this paper. That algorithm relied on a completeness relation and transformations 
of HO states with the help of HO brackets. Even though the algorithm of Ref. [36| required calculation of onc- 
dimcnsional radial integrals, while the the present algorithm requires evaluation of two-dimensional radial integrals, 
the present algorithm is substantially more efficient. 



III. CONVERGENCE TEST FOR :i H AND 4 HE 



In this section, we apply the matrix elements of the N 2 LO %EFT NNN interaction obtained in this paper to the 
NCSM calculation of 3 H ad 4 Hc ground state properties. As a test of correctness of the computer code, we verified 
that the new more efficient algorithm reproduces the results obtained using the algorithm of Ref. [36[ for the two- 
pion-exchange term matrix elements. For the contact terms, we verified that in the limit of A — > 00, the matrix 
elements (|15p and lead to the same result and the same is true for matrix elements ([55]) and ([50]) . In addition, 
we benchmarked the computer code for evaluation of (fl"4]) and ([30]) with the computer code written by A. Nogga [37l |. 
Finally, we tested numerically that the use of (|3"3"|) results in the same matrix element as the use of (|3"5|) in the three- 
nucleon antisymmetrized basis \NiJT) introduced in Eq. ||SJ). The same checks were also performed for the alternative 
version of the one-pion-exchange plus contact term (D-term) given in Eq. (|29p . That is, we verified numerically that 
the matrix element ([4T)|) leads to the same result as (f3"0")) in the limit of A — > 00 and the use of (|4"D|) results in the same 
matrix element in the three-nucleon antisymmetrized basis \NiJT) as the use of (|33|) multiplied by (— 1)* + * +s+s . 

We use the N 3 LO NN interaction of Ref. [lj|. We adopt the c\, C3 and C4 LECs values as well as the value of A 
from the N 3 LO NN interaction of Ref. [r| for our local chiral EFT N 2 LO NNN interaction. The regulator function 
was chosen in a form consistent with that used in Refs. [ll[ and [l3|: F(q 2 ; A) = exp(— q 4 /A 4 ) (fTTj) . Values of the Cd 
and ce LECs are constrained by a fit to the A = 3 system bind ing energ y fl4 . [l8|. Obviously, additional constraints 
are needed to uniquely determine values of cd and ce, see Refs. [ll|, [fl EM \2M for discussions of different possibilities. 
Here we are interested only in con verg ence properties of our calculations. Therefore, we simply select a reasonable 
value, e.g. cr> = 1, and follow Ref. [lj| and adopt ce value as an average of fits to 3 H and 3 He binding energies. In 
Table [II we summarize the NNN interaction parameters used in calculations described in this section. We note that 
4 He results obtained with the identical Hamiltonian but with cd = —1 and ce = —0.331 are presented in Ref. [l7| . 



TABLE I: NNN interaction parameters used in the present calculations. The regulator function was chosen in the form 
F(q 2 ;A)=c X p(-q 4 /A 4 ). 



ci [GeV"l] 


c 3 [GeV-1] 


c 4 [GeV"l] 


CD 


c E 


A [MeV] 


A x [MeV] 


AU[MeV} 


9A 


Ft, [MeV] 


-0.81 


-3.2 


5.4 


1.0 


-0.029 


500 


700 


138 


1.29 


92.4 



We use the Jacobi coordinate HO basis antisymmentrized according to the method described in Ref. [35[. In Figs.[(5] 
and [3 we show the convergence of the 3 H ground-state energy and point-proton rms radius, respectively, with the 
size of the basis. Thin lines correspond to results obtained with the NN interaction only. Thick lines correspond to 
calculations that also include the NNN interaction. The full lines correspond to calculations with two-body effective 
interaction derived from the chiral EFT N 3 LO NN interaction. The dashed lines correspond to calculations with the 
bare chiral EFT N 3 LO NN interaction. The bare NNN interaction is added to cither the bare NN or to the effective 
NN interaction in calculations depicted by thick lines. We observe that the convergence is faster when the two-body 
effective interaction is used. However, starting at about iV max = 24 the convergence is reached also in calculations 
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with the bare NN interaction. The rate of convergence also depends on the choice of the HO frequency. In general, it 
is always advantageous to use the effective interaction in order to improve the convergence rate. The 3 H ground-state 
energy and point-proton radius results are summarized in Table [TTJ The contributions of different NNN terms to the 
3 H ground-state energy are presented in Table ITTT1 In addition to results obtained using the cd = 1, we also show 
in Table InTl results obtained using cd = —1 and a corresponding eg constrained by the avarage of the 3 H and 3 He 
binding energy fit. For completeness, we show results obtained by the two alternative one-pion-exchange plus contact 
terms ([28)1 and (|29| . In all cases, the contact -E-term gives a positive contribution. The contribution from the D-teim 
changes sign depending on the choice of cz>. Still, the two-pion exchange c-terms dominate the NNN expectation 
value. 




A -A NN 28 bare 
>->NN28 eff 
O- O NN+NNN 28 bare 
<^> NN+NNN 28 eff 



2 



8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 
N 



FIG. 6: (Color online) 3 H ground-state energy dependence on the size of the basis. The HO frequency of HQ — 28 MeV was 
employed. Results with (thick lines) and without (thin lines) the NNN interaction are shown. The full lines correspond to 
calculations with two-body effective interaction derived from the chiral NN interaction, the dashed lines to calculations with 
the bare chiral NN interaction. For further details see the text. 



In Figs. [8] and [9l we show convergence of the 4 He ground-state energy and point-proton rms radius, respectively. 
The NCSM calculations are perforemed in basis spaces up to N max = 20. Thin lines correspond to results obtained 
with the NN interaction only, while thick lines correspond to calculations that also include the NNN interaction. The 
dashed lines correspond to results obtained with bare interactions. The full lines correspond to results obtained using 
three-body effective interaction (the NCSM three-body cluster approximation). It is apparent that the use of the 
three-body effective interaction improves the convergence rate dramatically. We can see that at about iV max = 18 the 
bare interaction calculation reaches convergence as well. It should be noted, however, that p-shell calculations with 
the NNN interactions are presently feasible in model spaces up to iV max = 6 or N max = 8. The use of the three-body 
effective interaction is then essential in the p-shell calculations. 

We note that NCSM calculations in the three-body cluster approximation are rather involved. The 4 He NCSM 
calculation with the three-body effective interaction proceeds in three steps. First, we diagonalize the Hamiltonian 
with and without the NNN interaction in a three-nucleon basis for all relevant three-body channels. In the second 
step, we use the three-body solutions from the first step to derive three-body effective interactions with and without 
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2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 

N 



max 

FIG. 7: (Color online) 3 H point-proton rms radius dependence on the size of the basis. The HO frequency of fr£l — 28 MeV was 
employed. Results with (thick line) and without (thin line) the NNN interaction are shown. The two-body effective interaction 
derived from the chiral NN interaction was used in the calculation. For further details see the text. 



TABLE II: Ground-state energy and point-proton rms radius of 3 H and 4 He calculated using the chiral N 3 LO NN potential 
[l^ | with and without the local chiral N 2 LO NNN interaction. The LECs values and other parameters are given in Table [T] 
The calculations were performed within the ab initio NCSM. 





NN 


3 H 


NN+NNN 


Expt. 


E gB [MeV] 


-7.852(5) 




-8.473(5) 


-8.482 


r p [fm] 


1.650(5) 




1.608(5) 








4 He 








NN 




NN+NNN 


Expt. 


E ss [MeV] 
r p [fm] 


-25.39(1) 
1.515(2) 




-28.34(2) 
1.475(2) 


-28.296 
1.455(7) 



TABLE III: Contributions of different NNN terms to the 3 H ground-state energy. The cd and ce LECs are explicitly shown. 
Other parameters are given in Table [I] All energies are given in MeV. The two alternative one-pion-exchange plus contact 
terms ([28} and (I29|) are considered. 





3 H 










CD 


ce 


E ga 


c terms 


D term 


E term 


1.0 (Eq. (28}) 


-0.029 


-8.473 


-1.01 


0.13 


0.03 


-1.0 (Eq. ([28}) 


-0.331 


-8.474 


-1.07 


-0.16 


0.32 


1.0 (Eq. (HI) 


-0.159 


-8.471 


-0.99 


0.005 


0.14 


-1.0 (Eq. ([29])) 


-0.213 


-8.474 


-1.10 


-0.05 


0.21 
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A-ANN 36 bare 

NN 36 eff 
l>-t>NN28 eff 
O- O NN+NNN 36 bare 
□— □ NN+NNN 36 eff 
<^> NN+NNN 28 eff 




FIG. 8: (Color online) 4 He ground-state energy dependence on the size of the basis. The HO frequencies of HQ, = 28 and 
36 MeV was employed. Results with (thick lines) and without (thin lines) the NNN interaction are shown. The full lines 
correspond to calculations with three-body effective interaction, the dashed lines to calculations with the bare interaction. For 
further details see the text. 



the NNN interaction. By subtracting the two effective interactions we isolate the NN and NNN contributions. This is 
needed due to a different scaling with particle number of the two- and the three-body interactions. The 4 He efffective 
interaction is then obtained by adding the two contributions with the appropriate scaling factors [T(|. In the third 
step, we diagonalize the resulting Hamiltonian in the antisymmetrized four-nucleon Jacobi-coordinate HO basis to 
obtain the 4 He J n T = + ground state. Obviously, in calculations without the NNN interaction, the above three 
steps are simplified as no NNN contribution needs to be isolated. In addition, in the case of no NNN interaction, 
we may use just the two-body effective interaction (two-body cluster approximation), which is much simpler. The 
convergence is slower, however, see discussion in Ref. f39[. We also note that 4 Hc properties with the chiral N 3 LO 
NN interaction that we employ here were calculated using two-body cluster approximation in Ref. [40| and present 
results are in agreement with results found there. 

Our 4 He results are summarized in Table |TTJ We note that the present NCSM 3 H and 4 He results obtained with 
the chiral N 3 LO NN interaction are in a perfect agreement with results obtained using the variational calculations 
in the hyperspherical harmonics basis as well as with the Faddeev-Yakubovsky calculations published in Ref. [4l| . 
A satisfying feature of the present NCSM calculation is the fact that the rate of convergence is not affected in any 
significant way by inclusion of the NNN interaction. 



IV. CONCLUSIONS 



In this paper, we regulated the NNN interaction derived within the chiral effective field theory at the N 2 LO with 
a function depending on the magnitude of the momentum transfer. The regulated NNN interaction is local in the 
coordinate space. This is advantages for some many-body techniques. In addition, it was found that this interaction 
performs sligthtly better in mid-p-shell nuclei than its nonlocal counterpart [HI, [37j . We calculated matrix elements 
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2 4 6 8 10 12 14 16 18 

N 



max 

FIG. 9: (Color online) 4 He point-proton rms radius dependence on the size of the basis. The HO frequencies of Ml = 28 and 36 
MeV was employed. Results with (thick line) and without (thin line) the NNN interaction are shown. The three-body effective 
interaction was used in the calculation. For further details see the text. 

of the local chiral NNN interaction in the three-nuclcon HO basis and performed calculations for 3 H and 4 Hc within 
the ab initio NCSM. We demonstrated that a very good convergence of the ground-state properties of these nuclei 
remains unchanged when the NNN interaction is added to the Hamiltonian. Expressions for the local x^FT NNN 
interaction matrix elements derived in this paper may by used after some modifications with other bases, e.g. with 
the hypersphcrical harmonics basis. 
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